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Abstract. - We outline a generic ratchet mechanism for creating directed spin-polarized currents 
in ac-driven double well or double dot structures by employing resonant spin transfer through the 
system engineered by local external magnetic fields. We show its applicability to semiconductor 
nanostructures by considering coherent transport through two coupled lateral quantum dots, where 
the energy levels of the two dots exhibit opposite Zeeman spin splitting. We perform numerical 
quantum mechanical calculations for the I-V characteristics of this system in the nonlinear regime, 
which requires a self-consistent treatment of the charge redistribution due to the applied finite 
bias. We show that this setting enables nonzero averaged net spin currents in the absence of net 
charge transport. 
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Introduction. — The field of semiconductor spintron- 
ips has seen rapid progress lately, yet there are still many 
obstacles on the way from fundamental research to oper- 
ating spin-based devices [1]. The creation of spin polar- 
ized currents is one basic requirement for the realization of 
semiconductor spintronics systems that share the prospect 
of being able to outperform conventional electronics. Due 
to a better controlability and faster processing times it is 
favorable to generate those currents by electrical means, 
e.g. by the variation of (contact) voltages. Promising 
classes of devices include spin pumps [2-5] , spin rectifica- 
tion [6] and spin ratchets [7-12]. These proposals share 
the common idea to generate directed spin currents, e.g. 
mediated by spin-orbit interaction, upon time variation 
of external potentials. Here we focus on spin ratchets, 
a generalization of the particle quantum ratchet mecha- 
nism [13-15]. In such systems with broken spatial sym- 
metry, pure spin currents are generated by means of an 
ac-driving with no net average bias. This idea has been 
put forward for both, nonlinear ly driven coherent conduc- 
tors [7-9], as well as conductors in the dissipative regime, 
where Brownian particle motion is converted into directed 
spin currents [10-12]. While a net spin current could be 
shown to exist for the different settings, its magnitude is 
difficult to predict and an optimization towards larger spin 
currents is often not evident. 

Here we propose another, generic, spin ratchet mecha- 
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nism that is based on coherent resonant charge and spin 
transfer. It thereby leads to larger and controllable output 
and can be implemented in a variety of systems. More- 
over, since the ratchet spin currents require operation un- 
der nonequilibrium conditions and since the spin currents 
can usually be enhanced for strong ac-bias, we employ 
a fully self-consistent treatment of the electrostatics for 
our quantum transport calculations in the nonequilibrium 
regime. This involves a self-consistent determination of 
the voltage drop across the ratchet, which has been ap- 
proximated so far only by simple heuristic models [7-9]. 

After outlining the general working principle we focus 
on a setup invoking resonant tunneling through two quan- 
tum dots (QD) in a two-dimensional electron gas (2DEG). 
In the literature double QD systems have already been 
proposed as spin filters [16, 17] or sources for pure spin 
currents [18]. However, these proposals are based on QDs 
in the Coulomb blockade regime, whereas the double QDs 
considered here are strong coupled to the leads, i.e. trans- 
port is fully coherent and the conductance is larger. 

Mechanism. — To illustrate the envisioned spin 
ratchet mechanism, let us first consider the simplified one- 
dimensional potential model shown in Fig. la). Three 
electrostatic barriers divide the system into four regions 
(R1-R4). While the regions Rl and R4 support states 
with a continuous energy spectrum, the regions R2 and 
R3, representing the double QD structure, accommodate 
discrete resonant states due to the confinement imposed 
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Fig. 1: (Color online) Illustration of the working principle of a resonant spin ratchet. Panel a): system devided into four regions 
Rl to R4 through three barriers. A magnetic field, oriented in opposite directions in regions R2 and R3, splits the levels of 
spin-up and spin-down electrons (one representative level shown). Panels b,c): Upon application of a positive or negative bias 
voltage, ±Vb, the transmission probability is resonantly enhanced for spin- up or spin-down electrons, respectively. 



by the barriers. We further assume a magnetic field ori- 
ented in opposite directions in R2 and R3. Thereby, the 
resonant energy levels in R2 and R3 are spin split due to 
the Zeeman coupling, however oppositely in both regions. 
A finite bias voltage across the device shifts the energy lev- 
els in R2 with respect to those in R3. This enables one to 
bring energy levels of the same spin state in the two QDs 
into resonance at a specific bias voltage, see Fig. lb,c). 
For charge transport through the device, this results in 
an enhanced transmission of electrons of that specific spin 
state. Considering both forward and backward bias, it is 
obvious that the energy levels of different spin states can 
be brought into resonance for different signs of the bias 
voltage (see Fig. lb,c)). Therefore, upon applying an ac- 
bias to the system, spin can be transported in the absence 
of a net charge current, as we will confirm below. 
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Fig. 2: (Color online) a) Possible experimental realization of 
the principle mechanism for spin current generation described 
in Fig. 1: Two quantum dots are electrostatically defined (e.g. 
via side gates) within a quantum wire patterned on a 2DEG. 
The magnetic fringe fields of two oppositely magnetized ferro- 
magnetic stripes create the opposite Zeeman splitting in the 
two quantum dots, b) Sketch of an additional electrostatic 
potential, see text. 



System and method. — In Fig. 2 we depict a possi- 
ble experimental setup, which takes advantage of the prin- 
ciple just described and can be realized with present day 
material processing techniques. It is based on a quantum 
wire (QW) patterned on a 2DEG in the (x, y)-plane, see 
Fig. 2a. Within this wire (in x-direction) which is con- 
nected to two non magnetic leads, two QDs possessing 
discrete energy levels are formed, e.g. via side gates. To 
realize the opposite Zeeman splitting inside these QDs, 
two ferromagnetic stripes (FMS) with opposite in-plane 
magnetization (M = ±My) are patterned on top of the 
semiconductor heterostructure. The fringe fields of the 
FMS give rise to a non-uniform magnetic field B(x,y) in 
the plane of the 2DEG. We note that the proposed setup 
is just one possible realization of the mechanism to gen- 
erate pure spin currents outlined in Fig. 1. Alternatively, 
one could for instance think of charge transport through 
resonant tunneling diodes, where the Zeeman splitting can 
be introduced, e.g., by layers of dilute magnetic semicon- 
ductors [19,20]. 

The quantum dynamics of electrons in the conductor is 



described by the single-particle Hamiltonian 

n = + 9 -^B( X , y).ff+ V(x, y), (1) 

where m* is the effective mass, g* the effective gyroscopic 
factor, fiB the Bohr magneton, and a is the vector of the 
Pauli spin operators. Orbital effects due to the magnetic 
field are accounted for by the vector potential A(x,y), 
which enters the momenta IL, (x, y) — pi~ eAi(x, y), while 
the Zeeman term hg* [1bB(x, y) ■ a couples the spin degree 
of freedom to the external magnetic field B(x,y) due to 
the FMS. For this setup B(x,y) can be evaluated - along 
with A(x,y) - using standard magnetostatics [21]. The 
potential V(x,y) includes the lateral confining potential 
which forms the QW and QDs, a possible potential offset 
due to an additional gate voltage (see Fig. 2b), as well as 
the electrostatic potential due to an applied (driving) bias 
voltage Vo between the left and right contact to which the 
QW is connected. The driving is assumed to be adiabatic. 

We neglect inelastic processes and assume phase coher- 
ent electron transport. The bias eVo = (J,l — MR induces 
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an electrical current, which we evaluate in the right lead. 
Specifically, the current of electrons with spin polarization 
a = ± (with respect to a quantization axis in y direction) 
can be written as 



dE Af(E;V )T a (E;V ) . 



(2) 



E C 



Here Eq denotes the energy of the conduction band edge, 



f(E,E F -e% )] is 



and Af(E-Vo) = [f(E,E F + e%) 
the difference between the Fermi functions in the two 
leads. In Eq. (2), the quantum transmission probability 
for electrons with spin a is given by 



T a (E;V ) 



N L JV R 

EE E 



\tncr.n' a >(E;V ) 



(3) 



where t nn y a ' is the amplitude for transmission from the 
scattering state (nV) in the left lead into the scattering 
state (na) in the right lead, with the summations running 
over the -/V L / R open transversal channels of the left /right 
lead. These amplitudes are evaluated by projecting the 
Green's function of the open system onto an appropriate 
set of asymptotic spinors defining incoming and outgoing 
channels. Making use of a real-space discretization of the 
Schrodinger equation [22] , the calculation of the S- matrix 
elements was made feasible by the implementation of a 
recursive algorithm for the calculation of Green's functions 
for spin-dependent transport [23]. 

Self-consistent numerical procedure. Since 
ratchet (spin) currents are expected for ac-driving with 
external bias voltages Vo in the nonlinear regime, the pro- 
file of the electrostatic potential, dropping across the de- 
vice, may play an important role for the working princi- 
ple of the ratchet device, and hence its treatment needs 
special care. In the following we outline our approach to 
the non-equilibrium quantum transport problem including 
the self-consistent determination of the electrostatic po- 
tential drop arising from the charge rearrangement in the 
nonlinear bias regime. We consider the classical electro- 
static potential V es (f) described by the Poisson equation, 
V 2 V cs (r) = ep(f)/(se ), with e = 15.15 for InAs. The 
charge density p(f) = —e[n(r) — n^r)] consists of both 
the density distributions of the electrons, n(f), and the 
donors, nd{r). The latter is usually not known a priori, 
while the electron density can be calculated from the lesser 
Green function Q < , see e.g. [22]: 



n (f) = ~^ J dEQ<(r,r;E). 



(4) 



In equilibrium the electrostatic potential is typically in- 
cluded in the effective confinement potential V = V con f + 
V® s which we modelled by a hard- wall potential [24]. 
Therefore, the equilibrium electrostatic potential, gov- 
erned by V 2 V c °(f) = -e 2 [n (r) — rid(r)] / (eeo) , with no(f) 
the electron density for zero-bias, need not be considered 



explicitly Thus, in the noncquilibrium situation with a fi- 
nite source-drain bias, only the change of the electrostatic 
potential SV CS = V cs — V® s due to the charge rearrangement 
5n = n — no of the electrons has to be included, similar as 
described in Ref. [25] . Hence the Poisson equation 



V 2 SV cs (r) 



e 

ee 



i(r) - no(r) 



(5) 



has to be solved. 

To account for the influence of the leads it is conve- 
nient to write the electrostatic potential as a sum, SV CS = 
Vbias + Ki [25] . Here, Vbias is the potential induced by the 
contacts ignoring the charges in the device. It is the solu- 
tion of the Laplace equation V 2 Vbias (^) = with bound- 
ary conditions Vbi as (x = ±L/2) = ±eVb/2, where L is 
the distance between the contacts. For three-dimensional 
contacts, Vbias has the shape of a linear ramp. The con- 
tribution of the charges inside the system is described by 
Vol, which is the solution of Eq. (5) with the boundary 
conditions of V c \(x = ±L/2) = at the interfaces of the 
contacts. 

If the electron density in the leads is substantially higher 
than in the device, the potential drop is screened in the 
leads leading to a flat electrostatic potential close to the 
contacts. Then, n(r') no(f ') far from the device, and 
we can compute the electrostatic potential as 



Aneeo 



2 -no(r') 



d z r 



\r — r' 



(6) 



which is the solution of the Poisson equation for a vanish- 
ing potential at infinity. 

For the self-consistent solution of the transport prob- 
lem we compute from Eq. (4)the electron density n{r) for 
a given electrostatic potential and obtain an improved 
profile for the electrostatic potential from n(r). It is 
known, however, that the straightforward iteration be- 
tween Eqs. (4) and (6) typically does not converge [26]. 
Instead we adapted the Newton-Raphson method intro- 
duced in Refs. [26, 27] in order to evaluate the electro- 
static potential by means of Eq. (6), which significantly 
improves the convergence behavior of the self-consistent 
computation scheme. After calculating the self-consistent 
electrostatic potential for each value of the source-drain 
voltage we are able to determine the transport properties 
of the considered ratchet device using Eq. (2). 

Numerical results. — In the inset of Fig. 3 we show 
the transmission T(E; 0) in the linear-response regime, 
Vo — > 0, at zero magnetic field for the device depicted 
in Fig. 2. The system's geometry parameters (specified in 
detail below) and energies are chosen such that the leads 
carry 3 to 4 transversal modes, while the point contacts 
only allow for (resonant) tunneling. Thus, we find sev- 
eral sharp transmission peaks corresponding to the dis- 
crete resonant energy levels of the two QDs. The peaks 
appear as doublets due to the inter-dot tunnel splitting. 
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Fig. 3: (Color online) Spin resolved currents I a for spin-up 
(a — +, black curve) and spin-down (a = — , red curve) 
as a function of applied bias voltage Vb for an InAs-based 
double dot device (Fig. 2) at Fermi energy Ef = 4.66 meV. 
The solid (dashed) lines denote the calculation for the self- 
consistent (model) electrostatic potential drop. Inset: Trans- 
mission T(E, Vb = 0) through the same system at zero mag- 
netic field in linear response. 

Upon applying a voltage we expect an asymmetric I-V 
characteristic for the spin resolved currents for Fermi en- 
ergies close to those resonant energy levels. To conduct ex- 
plicit calculations for the spin current of the device shown 
in Fig. 2 we have to fix several parameters. Taking InAs 
as the material where the 2DEG is built from, we have 
m* = 0.024m and g* = 15. The width of the QW is 
chosen to be W — 200nm, while at the point contacts, 
which are separated by 450nm, the width narrows down 
to 70nm. The two FMS of identical size have dimensions 
xq = 250nm, tjq = 800nm, zq = 250nm and are centered 
on top of the two point contacts at a distance 50nm above 
the 2DEG. For the magnetization of the stripes we chose 
fx M — 2.5T, a value well in reach using e.g. Dyspro- 
sium [28] . To allow for a high electron density in the leads 
in order to screen the electrostatic potential drop and to 
achieve a flat electrostatic potential in Rl and R4, we in- 
clude an additional confinement potential Uo — 2.92meV 
as shown in Fig. 2b) . 

The spin-resolved I-V characteristics which are ob- 
tained with the described procedure for the Fermi energy 
at the second resonant energy level, E-p = 4.66 meV (see 
inset of Fig. 3), are shown in Fig. 3 as solid lines (black 
and red line for spin- up and down, respectively). There we 
find that the device indeed acts as a spin selective current 
rectifier with dominant (> 90%) spin- up polarization for 
positive and spin-down polarization for negative bias volt- 
ages. This implies that when applying an adiabatically 
varying unbiased ac-voltage to the contacts the system 
only transports spins, yet no net charge, since the total 
charge currents average to zero. The pronounced maxi- 
mum (minimum) at ±Va reflects level alignment, i.e. the 
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Fig. 4: (Color online) Full lines: Self-consistent electrostatic 
potential drop V ea (x) for the two quantum dots bounded by 
points contacts at x — xi,X2,X3 and averaged over the trans- 
verse direction. It is shown for bias voltages Va/b as indicated 
in Fig. 3. The dashed line corresponds to the heuristic volt- 
age drop model (underlying the dashed spin current curves in 
Fig. 3). 



resonance condition sketched in Fig. 1. For comparison 
the dashed lines in Fig. 3 represent the corresponding I a - 
V curves arising from the aforementioned heuristic voltage 
drop model. It assumes a local voltage drop (Vb/3) at each 
of the three point contacts, i.e. the regions of maximum re- 
sistance. We find the same qualitative overall behavior for 
the spin currents, while the maximum of the resonant tun- 
neling current is reached at a higher value, Vb = Vb with 
eVe being approximately three times the Zeeman split- 
ting of the energy levels of the QD, as expected for this 
particular model. We performed further calculations (not 
presented here) which show that the degree of spin polar- 
ization, i.e. the efficiency of the device, increases with the 
Zeeman splitting in regions R2 and R3, due to the stronger 
separation of the spin — er energy levels at resonance for 
spin er. Accordingly, the optimal working condition, i.e. 
levels for one spin species at resonance, is then reached at 
higher voltages. 

Since the amount of voltage that drops at the central 
point contact determines the alignment of the energy levels 
of the two QDs, in Fig. 4 we study the profile of the self- 
consistent electrostatic potential, V es (x), along the trans- 
port direction for applied bias Vb = Va. (full black line) 
and Vb (full red line) with Va/b marked in Fig. 3. We see 
that due to the point contacts a step structure emerges, 
with its exact form depending on Vb and Ep. Compared 
to the heuristic voltage drop model (dashed line), we find 
a stronger voltage drop at the central point contact. This 
explains that in Fig. 3 the current maximum is reached 
at lower bias voltages and the currents are significantly 
smaller at high bias for the self-consistent current calcu- 
lation. Still, in view of Figs. 3 and 4, we can conclude 
that the heuristic voltage drop model constitutes a fair 
approximation to 14s for the system considered. 

Summary. — In conclusion we have presented a 
generic mechanism to produce pure spin currents in nanos- 
tructures by applying an ac electrical bias. It is based 
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on resonant transfer of spin-up and -down electrons in 
opposite directions and gives rise to pure ratchet type 
spin currents upon driving. Our results are based on self- 
consistent Keldysh Greens function transport in order to 
adequately describe the nonequilibrium condictions under 
which the system works. 
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